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ABSTRACT 

With the availability of active weather radar observations from space from the 
Precipitation Radar (PR) on board the Tropical Rainfall Measuring Mission (TRMM) 
satellite, numerous studies have been performed comparing PR reflectivity and derived 
rain rates to similar observations from ground-based weather radars (GR). These studies 
have used a variety of algorithms to compute matching PR and GR volumes for 
comparison. Most studies have used a fixed 3-dimensional Cartesian grid centered on the 
ground radar, onto which the PR and GR data are interpolated using a proprietary 
approach and/or commonly available GR analysis software (e.g., SPRINT, REORDER). 
Other studies have focused on the intersection of the PR and GR viewing geometries 
either explicitly or using a hybrid of the fixed grid and PR/GR common fields of view. 

For the Dual-Frequency Precipitation Radar (DPR) of the upcoming Global Precipitation 
Measurement (GPM) mission, a prototype DPR/GR comparison algorithm based on 
similar TRMM PR data has been developed that defines the common volumes in terms of 
the geometric intersection of PR and GR rays, where smoothing of the PR and GR data 
are minimized and no interpolation is performed. The PR and GR volume-averaged 
reflectivity values of each sample volume are accompanied by descriptive metadata, for 
attributes including the variability and maximum of the reflectivity within the sample 
volume, and the fraction of range gates in the sample average having reflectivity values 
above an adjustable detection threshold (typically taken to be 18 dBZ for the PR). 
Sample volumes are further characterized by rain type (Stratiform or Convective), 
proximity' to the melting layer, underlying surface (land/water/mixed), and the time 
difference between the PR and GR observ ations. 

The mean reflectivity differences between the PR and GR can differ between data sets 
produced by the different analysis methods; and for the GPM prototy pe, by the ty pe of 
constraints and categorization applied to the data. In this paper, we will show results 
comparing the 3-D gridded analysis "black box” approach to the GPM prototype volume- 
matching approach, using matching TRMM PR and WSR-88D ground radar data. The 
affects of applying data constraints and data categorizations on the volume-matched data 
to the results will be shown, and explanations of the differences in terms of data and 
analysis algorithm characteristics will be presented. Implications of the differences to the 
determination of PR/DPR calibration differences and use of ground radar data to evaluate 
the PR and DPR attenuation correction algorithms will be discussed. 
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1. INTRODUCTION 2. DATA AND ANALYSIS CHARACTERISTICS 


Numerous studies have compared reflectivity and 
derived rain rates from ihe space-based Precipitation 
Radar (PR) on board the Tropical Rainfall Measuring 
Mission (TRMM) satellite to similar observations from 
ground-based weather radars (GR). using a variety of 
algorithms to compute matching PR and GR volumes 
for companson. Most studies have used a fixed 3- 
dimensional grid centered on the ground radar (e g , 
Schumacher and Houze. 2000; Anagnostou et at.. 
2001; Liao et al., 2001; Wang and Wolff. 2009), onto 
which the PR and GR data are interpolated using a 
proprietary approach and/or commonly available GR 
analysis software (SPRINT, REORDER). Other 
studies have focused on the intersection of the PR 
and GR viewing geometries either explicitly (Boten 
and Chandrasekar. 2000), or using a hybrid of the 
fixed grid and PR/GR common fields of view. For the 
Dual-Frequency Precipitation Radar (DPR) of the 
upcoming Global Precipitation Measurement (GPM) 
mission, a protolype DPR/GR comparison algorithm 
based on TRMM PR data has been developed that 
defines the common volumes in terms of the 
geometric intersection of PR and GR rays, where 
smoothing of the PR and GR data are minimized and 
no Interpolation is performed (Schwaller and Morris. 
2011 ). 

The mean reflectivity differences between the PR and 
GR can differ between data sets produced by Ihe 
different volume matching methods: and for Ihe GPM 
prototype, by the type of constraints and 
categorization applied to the data. In Ihis paper, we 
will show results comparing the 3-D gridded analysis 
“black box” approach to the GPM prototype geometry- 
matching approach, using matching TRMM PR and 
WSR-88D ground radar data. The effects of applying 
data constraints and data categorizations on Ihe 
volume-matched data to the results, and explanations 
of the differences In terms of dafa and analysis 
algorithm characteristics are presented below. 
Implications of the differences to the determination of 
PR/DPR calibration differences and use of ground 
radar data to evaluate the PR and DPR attenuation 
correction algorithms are also discussed. 
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The geometry matching algorithm calculates PR and 
GR averages at the geometric intersection of the PR 
rays with the individual GR radar elevation sweeps. 
The along-ray PR data are averaged only In the 
vertical, between the top and bottom height of each 
GR elevation sweep i! mtersecls (Figure 1 ). GR range 
bins are horizontally averaged over an area of 
coverage defined by the half-power points of each PR 
ray intersected, distance-weighted from the parallax- 
adjusted center of Ihe PR beam. Each GR elevation 
sweep is treated separately. The volume-matched 
data are a set of conical surfaces retaining the vertical 
coverage defined by the elevation sweeps of the GR 
volume scan, but with horizontal resolution and 
location redefined by the PR's scan/ray coordinates. 
The data gaps between GR sweeps and the “cone of 
silence' above (he highest sweep angle are retained 
in Ihe geomelry-malch data set. 



Figure I. Schematic of PR ray GR sweep intersections. 
Shaded areas arc "matching volumes" showing the PR 
gates for one PR ray intersecting GR sweeps (dashed) at 
two different elevation angles. PR gates are 250 nt along- 
ray by -5 km in the horizontal. 

Unlike the gridded approaches there is no 
Interpolation, extrapolation, or oversampling of data. 
SO matching volumes only exist at somewhat random 
locations where both the PR and GR instruments 
have taken actual observations. However, other than 
lor Ihe averaging required to produce the matching 
volumes, the data are not smoothed; and each 
sample volume is accompanied by metadata 




describing the variability and maximum of the 
reflectivity within the sample volume, and the fraction 
of range gales In the PR and GR sample averages 
having reflectivity values above an adjustable 
detection threshold (typically taken to be 18 dBZ for 
the PR). Sample volumes are further characterized 
by rain type (Stratiform or Convective), proximity to 
the melting layer, underlying surface 
(land /water/mixed), and the time difference between 
the PR and GR observations. 

The approaches using analysis of PR and GR data to 
a fixed 3-dimensional grid centered on the GR treat 
the PR and GR data separately. While offering the 
simplicity of a regular coordinate system ol fixed 
location and size, grids represent the scan pattern of 
neither instrument and thus require some amount of 
smoothing, interpolation, and extrapolation to attempt 
to fill as many grid points as possible with data values 
and fill reasonable gaps in the GR volume scan. All 
resulting non-missing data points are treated equally, 
whether or not one or both instruments made 
observations In the volume represented by the grid 
box. 

In this study, we consider matched PR and GR 
reflectivity data from the grid-based volume matching 
algorithm and the geometry-match algorithm. PR 
data are from the TRMM 2A-25 attenuation-corrected 
reflectivity product. Version 6. GR data originate from 
the WSR-880 Level II Archive reflectivity product, 
which has been quality-conirolled to remove non- 
precipitating echoes (Wolff at al, 2005). Only data 
samples within 100 km of the ground radar and the 
overlap of the PR data swath are evaluated. The 3-0 
grids used are of 4-km horizontal resolution and 1.5- 
km vertical resolution, with 13 levels ceniered 
between 1.5 and 19.5 km height above the GR. 

PR data are analyzed to the grid following the 
melhods applied by Liao, et al (2001). Two different 
grid analysis methods are applied to the GR data. 
The first method takes the 2-km-resolution 2A-55 
standard TRMM GV product and reduces it to 4 km 
resolution, as in Liao, et al. (2001). The second 
method analyzes the Level-li data to the 4-km, 13- 
level grid using Ihe REORDER radar analysis 
software. For purposes of comparison to the gridded 
data, the geometry-match data are grouped into the 
same 13 vertical levels based on the midpoint height 
of each sample volume. A mean bright band height is 
computed for each coincident PR/GR rain case from 
information provided by the PR bright band detection 
algorithm, in order to subdivide the data by proximity 
to the bright band (above, within, or below). 

In computing Ihe mean reflectivity differences 
between the PR and GR. the matched volumes are 
subdivided into categories based on combinations ol 
the following attributes common lo both Ihe grid- 
based and geometry-match data sets: 


■ TRMM orbit number (defines date and lime of ihe 
event) 

■ GR site identifier 

■ height layer (13 layers. 1.5-19.5 km) 

■ proximity lo bright band: above, within, or below 

■ rain type: stratiform, convective, or unknown 

■ distance from the GR (0-50, 51-100 km) 

For each of the data categories defined by the 
permutations of these attributes, the mean difference 
between, and standard deviation of, Ihe PR and GR 
reflectivity for the non-missing sample volumes in Ihe 
category is computed separately for the grid data and 
Ihe geometry match data and stored in a data table, 
along with the identifying attributes and the number of 
data samples included in the category. 

Geometry-match data are subdivided by an additional 
attribute defined as the fraction of the sample with 
reflectivity above a minimum instrument detection 
threshold, defined as 18 dBZ for PR and 15 dBZ for 
the GR (to match the PR detection threshold but allow 
for a 3 dBZ calibration difference). The geometry 
matching algorithm determines, from a pure 
geometric standpoint. Ihe locations of the PR and GR 
range bins that are "coincident", and the number of 
each (number PR expected, number GV expected). 
Then the reflectivity values of each range gate are 
evaluated before averaging. The number of PR bins 
below Ihe 18 dBZ threshold (number PR rejected) and 
the number of GR bins below 15 dBZ (number GV 
rejected) are computed and related to each PR and 
GR sample volume To compute Ihe PR volume 
average, Ihe algorithm leaves out (hose range bins 
below 18 dBZ and averages the remaining (the same 
approach is taken in determining (he vertically- 
averaged PR reflectivity for a fixed layer in the grid- 
based algorithm). No range bins are left out in 
computing the reflectivity average, maximum, and 
standard deviation for the GR sample volumes, bul 
those bins below 0.0 dBZ are set lo 0.0 dBZ. 

From these attributes, a percentage o( each sample 
volume that is above its respective detection 
fhreshold is computed for the geometry-match PR 
and GR. Samples where both the PR and GR 
percenl-above-threshold is non-zero includes all data 
points with a non-missing reflectivity value, and is akin 
to the grid-based approach. Restricting the data to 
samples with a PR and GR percent-above-threshold 
constraint of 100% provides the best and fairest 
comparison between the PR and GR instruments, 
where the entire PR sample volume is above the PR 
detection threshold, and the entire GR sample volume 
is filled wilh echoes above the PR detection threshold. 
One of the major goals of this study is to show the 
effects ol varying the percent above threshold criteria 
on the PR-GR mean reflectivity differences. This 
sludy computed mean differences from the geometry- 
matched data for 1 1 categories of percent-above- 
Ihreshotd cutoff, ranging between 0 and 100%, by 
10% steps. 



3. SENSITIVITY TO FRACTION OF SAMPLE 
VOLUME ABOVE DETECTION THRESHOLD 

Figures 2-5 show mean PR-GR reflectivity differences 
for all rainy overpasses at the KMLB (Melbourne, 
Florida) WSR-88D site front 13 August 2006 to 30 
June 2008. KMLB was selected since previous 
studies have shown it to be closely calibrated lo the 
PR and to have a stable calibration over time (Liao et 
al„ 2001, Liao and Meneghmi, 2009a). Figure 2 
shows the differences lor the convective rain, above 
bright band category, where the differences based on 
the geometry-match data have been further 
subdivided on a sample-by-sample basis by their 
percent of gates above threshold as described in the 
previous section. Outside of the percent above 
threshold, the grid-based results are for the matching 
categories (orbits, slle, rain type, proximity to bright 
band) Data at height levels above the bright band 
are merged. Categories where no geometry match 
samples meet the percent above threshold criteria are 
eliminated from both the gridded and geometry match 
data for that percentage, but the gridded data are not 
otherwise filtered on a sample-by-sample basis. 


PR-GR Reflrctlviry Blnti, by Percent filled Threshold 
Convective | Above Bright tend | KMLB | 8/7006 6/1 DOB 
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will Increase wflh percent above threshold from just 
above 0 0 dBZ at percentage values above 0. to 15 
dBZ or greater a I for samples where 1 00% of the GR 
bins in Ihe average are above threshold. The mean 
differences computed from the gridded data takes all 
matched PR and GR grid points in the category where 
the reflectivity values for both are 1 8 dBZ or greater. 

Figure 2 shows that the PR is high biased relative lo 
the GR by about 2 dBZ in the grid-based analyses, 
and by 1 dBZ or less In the geometry-match analyses. 
The high bias of the PR relative to ihe GR in Ihe latter 
data lowers from 1 .26 dBZ lo 0. 16 dBZ in the 
geometry-match data as the percent-above-threshold 
constraint increases from 0 to 100 and the "floor" 
reflectivity for the GR sample volumes included 
increases lo 15 dBZ, closer to the PR cutoff at 18 
dBZ. The grid-based analyses do not change 
Significantly with the change in ihe percent threshold 
since the all sample volumes are included for each 
category. Minor changes occur where grid data for 
some orbits are excluded when the geometry-match 
data for the same orbit have no sample volumes 
meeting the percent-above-threshold criterion of the 
data category. 



Figure 2. PR-GR rellcciiviry mean differences ai KMLB 
for convective rain samples above the bright band, by 
percent above threshold category (sec text). PR-2A55 and 
PR-REORDER series arc based on gridded analyses. PR- 
GR series is from geometry-matched data, using percent 
above threshold categories from 0 lo 100%. 

Noie ihe effect of varying Ihe percent above threshold 
criteria on the PR and GR geometry-match results. 

As the percent of the sample volume filled with above- 
delection-fhreshold reflectivity bins increases, the 
high bias of Ihe PR to the GR decreases, and vice 
versa. Much of this Is explained by the averaging 
technique, where only PR bins ol 18 dBZ or greater 
are Included in the PR average, while for the GR. all 
bins are included in the volume average, though the 
GR percent above threshold measurement for the 
geometry-matched data is based on the fraction of the 
GR bins at 15 dBZ or greater. Thus, regardless of the 
percent above threshold critehon applied to the PR, 
the lowest PR reflectivity will always be 18 dBZ or 
grealer. The lowest possible geometry-match GR 
reflectivity included in the mean difference calculation 


Figure 3. As in Fig. 2, but for stratiform rain type. Order 
of series is changed from Fig. 2, lor visibility. 

Figure 3 shows the results for the stratiform rain, 
above bright band Category. In this case Ihe grid- 
based PR-GR bias based on Ihe 2A-55 GR product is 
smaller than the bias based on Ihe REORDER 
analysis of Lhe GR volume scan, while the geometry 
match data exhibits the same tendencies but slightly 
higher PR-GR bias than the convective case. The 
smaller mean reflectivity differences for Die grid- 
based results compared to the convective case are 
due lo the lower overall reflectivity in the stratiform 
rain areas, where imposition of an 18 dBZ minimum 
for the gridpoint sample volumes included in the mean 
difference calculation puts the grid data In situation 
approaching the 100% above-threshold constraint 
applied to the geometry- matched data. There is also 
likely to be some contamination of the bright band in 
the grid case, where the brighl-band-influenced data 
are filtered by excluding those fixed layers whose 
centers lie within 1000 m of the mean bright band. 



PR-GR Reflectivity Biases, by Percent Filled Threshold 
Stratiform | Above Bright Band | KMLB [ 8/2006-6/7 DC 4 



but, for greater ranges from the radar, the vertical 
extent of GR bins contributing to such layers may 
overlap the bright band, raising the GR reflectivity with 
respect to the PR, The actual top and bottom of each 
geometry-match sample volume is compared to the 
mean bright band height when determining whether 
the sample volume Is above, below, or affected by the 
bright band, so bright band contamination is less likely 
for these data. 

It Is this category (slraliform. above bright band) that 
Is used to evaluate calibration differences between 
the PR and ground radars, as attenuation of (he PR at 
Ku band is at its minimum, and strong horizontal 
gradients of reflectivity are no! present, minimizing the 
non-uniform beam filling effects. Figure 2 shows that 
the calibration offset is highly sensitive lo the method 
used to calculate matching PR and GR sample 
volumes, as well as to the parameters used to select 
the data samples included in the calculations. 


PR spreading of the high PR reflectivities over a wider 
area than they are observed, resulting in high PR 
reflectivities being differenced against lower GR 
reflectivities. Since the PR data are averaged only in 
the vertical in the geometry-malch analysis, this 
source of bias is not present in Ihese data. 

PR-OR R*fT«ct1vitv Bunt, by Pert#n* fclted Threshold 
Stratiform | B*iow 8;>fhl Band | KMlfi j BJ 2006-6/2068 



Figures 4 and 5 show the mean differences below the 
bright band for ihe convective and stratiform rain rate 
categories, respectively. The stratiform case In Fig. 5 
follows a similar Irend to Ihe above-bright-band 
categories with respecl to ihe change with percenl 
above threshold and the relative biases of (he Ihree 
data sets. The geometry match data for the 
convective case in Fig. 4 break Ihe patiern of 
monotonically decreasing PR-GR biases with 
Increasing peroenl above threshold. In this category, 
the PR and GR reflectivities change in a similar 
manner with percenl above threshold, perhaps due to 
ihe attenuation corrections applied to the PR data. 



PR-GR Reflectivity by Percent Filled Threthold 
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Figure 4. As in Fig. 2. but for convective samples below 
ihe bright band. 

The overall high bias of ihe convective samples for 
ihe gridded analyses relative to the geometry match 
data is due primarily to a few cases of very high 
convective reflectivities The mean reflectivity 
differences are weighted by Ihe number of gndpoints 
in the caiegory. not case-by-case, so a few cases with 
high PR-GR biases over targe areas are driving up 
the grid-based biases. The difference between the 
gridded data and geometry match dala in these cases 
is due to Ihe objective analysis scheme used for Ihe 


Figure S. As in Fig. 2, but for stratiform rain samples 
below the bright band. 

4. SENSITIVITY TO PR-GR TIME DIFFERENCES 

The lime matching rule for PR and GR data selects 
the GR volume scan with the earliest begin time in a 
9-minute window centered on the time of Ihe PR's 
closest approach to the GR site. The lime offset 
between the PR and GR dala has little effect on the 
mean reflectivity differences, as the mean PR and GR 
reflectivities do not change significantly in the range of 
lime offsets resulting from this rule. However, the 
point-to-point reflectivity differences for fast-moving or 
evolving precipitation echoes should be expected to 
increase as the time difference increases. To 
investigate these differences, the standard deviation 
of the point-to-point differences was computed for 
each caiegory, and averaged over Ihe full data set. 
Figure 6 shows Ihese results for Ihe gridded and 
geometry-malch data seis, for Ihe 100% above- 
threshold category. 



Figure 6. Standard Deviation of PR-KMLB reflectivity 
differences by time offset between PR and KMLB, for all 
categories shown in Figs. 2-5. combined. 


The standard deviation of the reflectivity differences 
increases for all 3 data sets as the time difference 




between the PR and GR increases from 0 to 3 
minutes. The reduction In the standard deviation at 4 
minutes time offset is probably a dala sampling effect 
due to the smaller number of data points in this group. 

5. SENSITIVITY TO MINIMUM REFLECTIVITY 
THRESHOLD 

By default, the geometry matching algorithm uses a 
PR threshold of 18 dBZ and a GR threshold of 15 dBZ 
in determining ihe fraction of a volume filled with 
above-threshold reflectivity. The sensitivity of the 
mean reflectivity differences to changes in these 
threshold values is demonstrated by changing Ihe GR 
threshold to 18 dBZ. to match the PR threshold. 
Table 1 shows mean PR-GR reflectivity differences 
for the two GR thresholds, split out into stratiform and 
conveclive rain regimes both above and below the 
bright band, limited to those samples t00% filled with 
above-threshold reflectivity. The data include all 
cases in years 2008 and 2009 ai KMLB As 
expected, the PR-GR mean differences for the 18 
dBZ GR threshold are lower than (or the 15 dBZ 
threshold, but only by about 0.3 (0.1) d8Z above 
(below) the bright band, and fewer samples (N) 
qualify (or the higher GR ihreshold. 

Tabic t. PR-KMLB mean reflectivity differences (dBZ) 
lor 200 & and 2009 from geometry-match data with GR 
reflectivity thresholds of 15 dBZ and 18 dBZ. Separate 
results are shown for convective (C) and stratiform (S) rain. 
above and below the bright band (BB). 


Rain Type 1 
Location 

15 dBZ GR 
Ihreshold 

16 dBZ GR 
threshold 

mean 

PR-GR 

N 

mean 

PR-GR 

N 

C / Above BB 

0.27 

1922 

-001 

1269 

C 1 Below BB 

1.03 

1154 

0.92 

1006 

S / Above BB 

-0.27 

2894 

-0.63 

1566 

S 1 Below BB 

2.17 

3174 

2.10 

2382 


6. SENSITIVITY TO RANGE FROM GR 

Table 2 shows the PR-GR geometry match mean 
reflectivity differences for KMLB for the dala periods 
used in Figs. 2-6, divided Into range categories of 0- 
50 and 50-100 Km from the GR. The sense In which 
the differences change with distance reverses 
between stratiform rain, where the differences 
increase with distance, and convective rain, where the 
differences decrease with increasing distance. The 
reason for this difference in behavior is nol 
immediately clear, as both the PR and GR volume 
averages are affected by the increase in the GR 
range gate height and width with distance. In either 
case, away from the bright band the difference 
between near and far distances in less lhan 0.4 dBZ 
for both convective and stratiform rain. The cause of 
the large differences with distance for the withm- 
bnght-band categories needs further investigation but 
may be a sampling issue due to (he smaller number 
of samples in the 0-50 km category. 


Table 2. PR-KMLB mean reflectivity differences (dBZ) 
for the geometry match data included in Figs. 2-6, split out 
by distance from the GR. Separate results are shown for 
conveclive (C) and stratiform IS) rain, above, below, and 
within the bright band (BB). 



0-50 km 

50-100 km 

Rain Type 1 
Location 

mean 

PR-GR 

N 

mean 

PR-GR 

N 

C 1 Above BB 

0 30 

165 

0 14 

1182 

C / Below BB 

1.55 

445 

1.17 

443 

C / Within B3 

3.04 

85 

0.37 

840 

S / Above BB 

-0 03 

237 

0.28 

1497 

S / Below BB 

1 19 

1540 

1.53 

1100 

S / Within BB 

-240 

105 

-0.66 

2818 


7. EFFECTS OF S-Ku FREQUENCY MATCHING 
ADJUSTMENTS 

All the comparisons shown up to this point have 
matched Ku-band PR reflectivity against S-band GR 
reflectivity, not accounting for expecled reflectivity 
differences due to Ihe different operating frequencies 
of each instrument. Liao and Meneghini (2009b) 
provide S- to Ku-band reflectivity corrections for the 
ice phase (above bright band) and rain phase (below 
bright band) based on theoretical considerations. 
Table 3 shows the results obtained comparing the 
geometry-match unadjusted (S-band) and Ku- 
adjusted GR reflectivities against the PR, for Ihe 
same data period as in Table 2 and Figs. 2-6. Note 
that no correction is attempted for the wilhin-bright- 
band layer, due to Ihe unknown particle sizes and 
types in this layer. 

Table 3. PR-KMLB mean reflectivity differences (dBZ) 
fur the geometry match data in Tabic 2. for both unadjusted 
and frequency-adjusted GR. Separate results arc shown for 
convective (C) and stratiform (S) rain, above and below die 


Rain Type / 
Location 

Unadjusted GR I Ku-adjusted GR 

mean 

PR-GR 

N 

mean 

PR-GR 

N 

C / Above BB 

0.16 

1347 

1.35 

1347 

C / Below BB 

1.36 

888 

-0.30 

888 

S / Above BB 

0.24 

1734 

0.73 

1734 

S / Below BB 

1 33 

2640 

0.61 

2640 


Note that the stratiform rain areas bolh above and 
below the melting layer show almost identical PR-GR 
mean reflectivity differences after the S-to-Ku GR 
adjustment. The S-to-Ku adjustment relationships are 
quadratic in lerms of Z v , Ihe reflectivity factor, 
resulting in larger adjustments to the convective 
cases. Assuming that the straliform/above bright 
band difference represents the residual calibration 
offset between the PR and GR. then applying this 
offset to the Ku-adjusted differences shows stratiform 
differences of 0.1 dBZ or less between PR and GR. A 
mean PR bias of approximately -1.0 dBZ exists for 
convective cases below the bright band, indicating an 
undercorrection for attenuation of the Version 6 PR at 
low levels in conveclive rain where PR attenuation is 




significant. These results are similar to those 
computed by Liao and Meneghini (2009b) lor KMLB, 
for post-orbital-boost cases between September 2001 
and February 2004. 

8. CASE-BY-CASE VARIABILITY 

Statistics shown thus far represent averages over all 
the cases in the time period. For comparison. Table 4 
presents mean PR-GR differences on a case-by-case 
basis (a raining TRMM overpass of the KMLB radar), 
for the stratiform rain, above bright band category, 
limited to those points with a percent above threshold 
of 100%. The results are ordered by the mean value 
of the maximum PR reflectivity in each remaining con- 
fixed sub-category (height and dislance in this case) 
and secondarily by orbit number. These data run 
from August 2006 to June 2008, as in Figs. 2-6. As 
seen in the results, the mean PR-GR differences for 
the geometry match data are insensitive to the mean 
reflectivity, with the exception of two outlier cases for 
orbits 60537 and 59408. However, the number of 
samples in the cases tends to increase with the 
maximum observed reflectivity in stratiform rain. 

Table 4. Casc-by-casc PR-K.MLB mean reflectivity 
differences fdBZ) for stratiform rain, above the bright band, 
PR-2AS5 and PR-RLORDER results are based on gridded 
PR and GR analyses. PR-GR Geo. Match results are from 
geometry-matched data, lor the 100% above threshold 
category. 


Orbit it 

Mean 

Max. 

PR 

PR-2A55 ' PR- 

g ridded REORDER 

gridded 

PR-GR 
Geo. Match 

Mean 

DIM. 

N 

Mean 

Dlff. 

N 

Mean 

Diff. 

N 


22 

0 76 

30 

KtMJ 

18 

-0 08 

5 


22 

1.28 

BFS1 


K1 

■cma 

6 

eeh 

23 

-0,48 

K3 


ma 


m 

cm 

msm 

0.16 

37 

■m 

■g 


8 


23 

0.48 

149 



-0.27 

EH 

50249 

24 




KETil 



56019 

24 


El 


S3 


El 

54691 

24 

0.33 

MB 


ma 

-0.41 

ma 

52676 

25 


ma 

-0.77 

mi 

1 44 

■Q 

50234 

25 

0.15 

msi 

-0.25 

ma 

-0.15 

ma 


25 


■a 

■syg 

ua 

BUI 

■u 


25 

0.34 

88 


. 76 

-0.50 
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Table 5 presents the case-by-case results for the 
convective rain, above bright band category. For this 
subset of data a pair of strong outlier cases appear for 
orbits 60537 and 51916 for all three analysis types. It 
is the large biases and numbers of samples for these 
cases that contribute to the high values of the PR-GR 
mean reflectivity differences for convective rain seen 
in the preceding figures. The reasons lor these outlier 
cases is a subject for further study. 


Table 5. As in Tabic 4, but for convective rain, above the 
bright band, 
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It is clear from Tables 4 and 5 that the case-by-case 
variability in the mean reflectivity difference between 
the PR and GR exceeds that of the effects of sample 
percent above threshold, minimum GR reflectivity 
threshold, range from the GR, and S-to-Ku frequency 
adjustments, and not all of this variation can be 
ascribed to the size of the data sample in each case. 

9. CONCLUSIONS 

A new volume-matching algorithm to compare space- 
based and ground-based radar observations has 
been developed for the upcoming GPM mission. Ii 
allows comparisons to be limited lo locations where 
boih sysiems observe echoes, with no Interpolation or 
extrapolation of the data, and allows the quality of the 
matching volumes to be controlled in terms of beam 
filling aspects. The geometry-matched data from this 
algorithm are compared lo traditional gnd-based 
analyses of the same data and are shown to produce 
a closer comparison between the TRMM PR and the 
Melbourne, Florida WSR-88D radar. 




































































































The two attributes that most affect the geometry- 
match comparison results are shown to be the 
percent of the matching volumes filled with reflectivity 
values above the PR detection threshold of 
approximately 18 dBZ. and the application o( S- to 
Ku-band frequency adjustments to the ground radar 
data, each of which can change the long-term mean 
reflectivity differences by up to 1.5 dBZ. Geometry- 
match and gnd-based comparison results for 
stratiform rain were similar, however for convective 
ratn the PR was much more high-biased against the 
GR lor gridded analysis when compared to the 
geometry-match result. 

Mean reflectivity differences were relatively 
insensitive to the lime difference between the PR and 
GR for the range of time differences allowed in the 
data set, though the scatter of the point-to-point 
differences is seen to increase with increasing lime 
differences. Mean PR-GR reflectivity differences as a 
function of distance from the ground radar trended in 
opposite directions for stratiform and convective rain, 
with a maximum absolute difference of about 0.4 d8Z 
for each. The case-by-case variability of the mean 
reflectivity differences was shown to exceed the 
variability In the full data set's differences resulting 
from any of the data analysis, categorization, and 
frequency adjustment methods applied In Ihe sludy. 

10. RESOURCES 

Time-malched TRMM PR and KMLB WSR-88D data 
files in original formats, geometry match netCOF data 
files produced from these dala, and the Data User's 
Guide for ihe geometry match data are freely 
available for download, as is open source code used 
to perform the geometry matching and generate 
displays and statistical comparisons between the PR 
and GR. Refer to ihe online links within the Validation 
Network Software and Data Products seclion of: 

http://pmm.nasa.gov/sciencefground-validation 
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